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We numerically study influence of a polychromatic perturbation on wave packet dynamics in one- 
dimensional double-well potential. It is found that time-dependence of the transition probability 
between the wells shows two kinds of the motion typically, coherent oscillation and irregular fluctu- 
ation combined to the derealization of the wave packet, depending on the perturbation parameters. 
The coherent motion changes the irregular one as the strength and/or the number of frequency com- 
ponents of the perturbation increases. We discuss a relation between our model and decoherence 
in comparing with the result under stochastic perturbation. Furthermore we compare the quantum 
fluctuation, tunneling in the quantum dynamics with ones in the semiquantal dynamics. 

PACS numbers: 05.45.-a, 03.65.-w 
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I. INTRODUCTION 

Quantum fluctuation can suppress chaotic motion of 
wave packet in the phase space due to the quantum in- 
terference, as seen in kicked rotor 0,0]. On the contrary, 
the quantum fluctuation can enhance the chaotic motion 
of wave packet due to tunneling effect as seen in kicked 
double- well model 0, 3 . The relation between chaotic 
behavior and tunneling phenomenon in classically chaotic 
systems is interestingand important subject in study of 
quantum physics [E Til- Recently, the semiclassical de- 
scription for the tunneling phenomena in a classically 
chaotic system have been developed by several groups 



Lin and Ballentine studied interplay between the tun- 
neling and classical chaos for a particle in a double-well 
potential with oscillatory driving force [13 . They found 
that coherent tunneling takes place between small iso- 
lated classical stable regions of phase space bounded by 
Kolmogorov-Arnold-Moser (KAM) surfaces, which are 
much smaller than the volume of a single potential well. 

Hanggi and the coworkers studied the chaos-suppressed 
tunneling in the driven double-well model in terms of the 
Floquet formalism [Tl|. They found a one-dimensional 
manifold in the parameter space, where the tunneling 
completely suppressed by the coherent driving. The 
time-scale for the tunneling between the wells diverges 
because of intersection of the ground state doublet of the 
quasienergies. 

While the mutual influence of quantum coherence and 
classical chaos has been under investigation since many 
years ago, the additional effects caused by coupling the 
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chaotic system to the other degrees of freedom (DOF) 
or an environment, namely decoherence and dissipation, 
have been studied only rarely 0, 0] as well as the 
tunneling phenomena in the chaotic system. Since mid- 
eighties there are some studies on environment-induced 
quantum decoherence by coupling the quantum system 
to a reservoir [Tj, [lj, [TJj . Recently quantum dissipation 
due to the interaction with chaotic DOF has been also 

studiedHHHIl!. 

In this paper we numerically investigate the relation 

quantum fluctuation, tunneling and decoherence com- 
bined to the derealization in wave packet dynamics in 
one-dimensional double-well system driven by polychro- 
matic external field. 

Before closing this section, we refer to a study on a 
delocalization phenomenon by a perturbation with some 
frequency components in the other model. Casati et al. 
have reported that the kicked rotator model with a fre- 
quency modulation amplitude of kick can be mapped 
to the tight-binding form (Loyld model) on higher- 
dimensional lattice in solid-state physics under very spe- 
cific condition [T^ . l2fH . Then the number M of the incom- 
mensurate frequencies corresponds the dimensionality of 
the tight-binding system. The problem can be efficiently 
reduced to a localization problem in M + 1 dimension. 
As seen in the case of kicked rotators, we can also expect 
that in the double-well system the coupling with oscilla- 
tory perturbation is roughly equivalent to an increase in 
effective degrees of freedom and a transition from a lo- 
calized wave packet to delocalized one is enhanced by the 
polychromatic perturbation. The concrete confirmation 
of the naive expectation is one of aims of this numerical 
work. 

We present the model in the next section. In Sect. 3, 
we show the details of the numerical results of the time- 
dependence of the transition probability between the 
wells based on the quantum dynamics. Section 4 contains 
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the summary and discussion. Furthermore, in appendix 
A, we gave details of the classical phase space portraits in 
the polychromatically perturbed double-well system and 
some considerations to the effect of polychromatic per- 
turbation. In appendix B, a simple explanation for the 
perturbed instanton tunneling picture is given. 



II. MODEL 

We consider a system described by the following 
Hamiltonian, 

H{t) = T + j-^Wy, (i) 
1 M 

A(t) = a- — Y^e t 31X1(^1 + 9 t ). (2) 

For the sake of simplicity, {e^} and {6i} are taken as 
6i = e, 8i — 0, i — 1, 2, .., M in the present paper. Then 
M is the number of frequency components of the exter- 
nal field and e is the perturbation strength respectively. 

are order of unity and mutually incommensurate 
frequencies. We choose off-resonant frequencies which 
are far from both classical and quantum resonance in the 
corresponding unperturbed problem. The parameter a 
adjusts the distance between the wells and we set a = 5 
to make some energy doublets below the potential bar- 
rier. Note that Lin et al. dealt with a double-well system 
driven by forced oscillator (Duffing- like model), there- 
fore, the asymmetry of the potential plays an important 
role in the chaotic behavior and tunneling transition be- 
tween the symmetry- related KAM tori [HI El- However, 
in our model the potential is remained symmetric during 
the time evolution process, and different mechanism from 
the forced oscillation makes the classical chaotic behavior 

00,111. 

In the previous paper [24| we presented numerical re- 
sults concerning a classical and quantum description of 
the field-induced barrier tunneling under the monochro- 
matic perturbation (M = 1). In the unperturbed double- 
well system (M — 0) the instanton describes the coherent 
tunneling motion of the initially localized wave packet. It 
is also shown that the monochromatic perturbation can 
breaks the coherent motion as the perturbation strength 
increases near the resonant frequency in the previous pa- 
per. In the classical dynamics of our model, outstanding 
feature different from previous studies is parametric in- 
stability caused by the polychromatic perturbation. 

Based on our criterion given below, we roughly esti- 
mate the type of the motion, i.e. the coherent and irreg- 
ular motions, in a regime of the parameter space spanned 
by the amplitude and the number of frequency compo- 
nents of the oscillatory driving force. It is suggested 
that the occurrence of the irregular motion is related to 
dissipative property which is organized in the quantum 
physics |25( . The classical phase space portraits and sim- 



ple explanation of relation to the dissipative property are 
given in appendix A. 



III. NUMERICAL RESULTS 

We use Gaussian wavepacket with zero momentum as 
the initial state, which is localized in the right well of the 
potential. 

iftq, t = 0) = (an) 1 ^ cxp{- (g ~J° )2 }, (3) 

where q$ = y/a ~ 2.236 means a bottom of the right 
well. The Gaussian wavepacket can be approximately 
generated by the linear combination of the ground state 
doublet as ip(q,t = 0) ~ -^((po(q) + (q) ) , where <po 
and if i denote the ground state doublet. The recurrence 
time for the wavepacket is T = 27rfi,/Aeoi ~ 9.4 x 10 3 
in the unperturbed case (M = 0), where Aei is the 
energy difference between the tunneling doublet of the 
ground state. We set the spread of the initial packet a = 
l/3.4(~ 0.3) and fi = 1.0 for simplicity throughout this 
paper. Indeed, the ammonia molecule is well described 
by two doublets below the barrier heigth in unperturbed 
case. 

We numerically calculate the solution ift(q,t) of time- 
dependent Schrodinger equation by using second order 
unitary integration with time step 5 ~ 10~ 2 . We define 
transition probability of finding the wave packet in the 
left well, 

P L (t)= [ mq,t)\ 2 dq. (4) 

J — oo 

In the cases that the perturbation strength is relatively 
small, Pl (t) can be interpreted as the tunneling probabil- 
ity that the initially localized wave packet goes through 
the central energy barrier and reaches the left well. We 
can expect that the transition probability is en- 

hanced as the number M of the frequency components 
increases up to some extent because of the increasing of 
the stochasticity in the total system. 

Figure 1 shows the time-dependence of Pl (t) for vari- 
ous combinations of e and M. Apparently we can observe 
the coherent and irregular motions. The coherent motion 
of the wave packet can be well-described by the semi- 
quantal picture in a sense that the wave packet does not 
delocalize to the fully delocalized state. The semiquantal 
picture decomposes the motion of the wave packet into 
evolution of the centroid motion and the spreading and 
squeezing of the packet Q. (See Subsect.3.5.) For ex- 
ample, in cases of relatively small perturbation strength 
(e = 0.4), coherent motion remains still up to relatively 
large M(= 5). 

It is important to emphasize that the tunneling contri- 
bution to the transition probability Pl is not so signifi- 
cant for large e and/or M . Then Pl may be interpreted 
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FIG. 1: (Color online) Tunneling probability Ph{i) as a func- 
tion of time t for various M's. (a)e = 0.4. (b)e = 0.8. The 
calculation time is same order to Heisenberg time in the un- 
perturbed case. 



as a barrier crossing probability due to the activation- 
transition because the energy of wave packet increases 
over the barrier height in the parameter range. Espe- 
cially, in the relatively large perturbation regime we can 
interpret the delocalized states as chaos-induced dereal- 
ization in a sense that the classical chaos enhances the 
quantum barrier crossing rate quite significantly. The 
chaotic behavior in the classical dynamics is given in ap- 
pendix A, based on the classical Poincare section and 
so on 26] . In the present section, we mainly focus on 
the transition of the quantum state from the localized 
wavepacket to delocalized state based on the data of nu- 
merical calculation. 



A. decoherence of the dynamics 

Once the wave packet incoherently spreads into the 
space as the M and/or e increase, the wavepacket is 
delocalized and never return to Gaussian shape again 
within the numerically accessible time. Apparently, we 
regard the delocalized quantum state as a decoherent 
state in a sense that the behavior of the wave packet 
is similar to that of the stochastically perturbed case. 



(See Fig. 5(a).) In case of relatively small perturbation 
strength (e = 0.4), the decoherence of quantum dynam- 
ics appears at around M ~ 7, and Pi(£) fluctuates ir- 
regularly in case of large M(= 10). In short, the irre- 
versible derealization of a Gaussian wave packet gen- 
erates a transition from coherent oscillation to irregular 
fluctuation of -Pi(i). We have confirmed that the sim- 
ilar behavior is also observed for other sets of values of 
the frequencies and the different initial phases of the 
perturbation. 

Here, we define a degree of coherence APl of the time- 
dependence of Pl (t) , based on the fluctuation of the tran- 
sition probability in order to estimate quantitatively the 
difference between coherent and incoherent motions. 

AP L ^^((P L (t)-(P L (t)) T ) 2 )T, (5) 

where (...)t represents time average value for a period 
T = 9.4 x 10 3 . Note that we used APl in order to express 
the decoherence of the tunneling osccilation of the transi- 
tion probability in the parametrically perturbed double- 
well system. On the other hand, the other quantities 
such as purity, linear entropy and fidelity, are sometimes 
used to characterize the decoherence of the quantum sys- 
tem [22) . The transition of the dynamical behavior based 
on the fidelity for description of the decoherence in the 
double- well system will be given elsewhere [2(|. 
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FIG. 2: Perturbation strength e dependence of the degree of 
coherence APl of tunneling probability for various M's. The 
APl is numerically estimated by Ph(t). 

Figure 2 shows the perturbation strength dependence 
of APl for various M's. We roughly divide the type of 
motion of wave packet into three ones as follows. In the 
coherent motions, the value of APl's is almost same to 
the unperturbed case, i.e. APl > 0.3, in which cases the 
instanton-like picture is valid |2g|. A simple explanation 
of the perturbed instanton is given in appendix B. 

In the irregular motions which are similar to the 
stochastically perturbed case, the value of AP^'s be- 
comes much smaller, i.e. APl < 0.2. As a matter of 
course, there are the intermediate cases between the co- 
herent and the irregular motions, 0.2 < AP^ < 0.3. Note 
that the exact criterion of the intermediate motion is not 
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important in the present paper because we can expect 
that the transitional cases approach to the irregular case 
in the long-time behavior. It should be stressed that 
the critical value e c exists, which divides the behavior of 
Ph(t) into regular and irregular motions. 
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FIG. 3: (Color online) Approximate phase diagram in the pa- 
rameter space spanned by the amplitude and the number of 
frequency components of the oscillatory driving force. The 
motions are plotted as different marks based on the degree 
of coherence APl- Circles(0)i crosses(x), and triangles(A) 
denote coherent motions (APl > 0.3), irregular motions 
(AP L < 0.2), and the transitional cases (0.2 < AP L < 0.3) 
respectively. 

Figure 3 shows a classification of the motion in the 
parameter space which is estimated by the value of the 
degree of coherence AP L . It seems that two kinds of the 
motion, i.e. coherent and irregular motions, are divided 
by the thin layer corresponding to the "transitional case" . 
As M increases, decoherence of the motion appears even 
for small e. The numerical estimation suggests that there 
are the critical values e c (M) of the perturbation strength 
depending on M. When the perturbation strength e ex- 
ceeds the critical value e c (M) for some M, the tunneling 
oscillation loses the coherence. The approximated phase 
diagram roughly same as the diagram generated by max- 
imal Lyapunov exponent of the classical dynamics. (See 
appendix A.) 



B. reduction of the tunneling period for e < e c 

In this subsection we give a consideration to the reduc- 
tion of the tunneling period in the regular motion regime 
e < e c . 

Figure 4 shows the e— dependence of the the period 
Ttun of the tunneling oscillation estimated by the numer- 
ical data Ph{t) hi the coherent motion regime e < e c . 
We can observe the monotonically decreasing of the tun- 
neling period as the perturbation strength increases. In 
the monochromatically perturbed case, the reduction of 
the tunneling period can be interpreted by applying the 
Floquet theorem to the quasi-energy states and the quasi- 
energy as the Hamiltonian is time-periodic H ~(t-\-2ir / T2) = 
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FIG. 4: Perturbation strength dependence of tunneling pe- 
riod Ttun for M = 1 and M = 2 in the coherent motion regime 
e < e c . 



H(t). When the wave packet does not effecively ab- 
sorb the energy from the external perturbation the time- 
dependence of the quantum state can be described by 
the linear combination of a doublet of quasi-degenerate 
ground states with opposite parity because we prepare 
the initial state in ip(q,t = 0) ~ -^75 (</?o(<?) + fiil)) an d 
the evolution is adiabatic. In the two-state approxima- 
tion that the avoided crossing of the eigenvalues dynam- 
ics does not appear during the time evolution, it is ex- 
pected that the state evolves as, 

TP(q, t)*± (u (q, t)e-> K E ^ dT + Ul (q, t)e~ l K £l ^ dT ) ,(6) 

where (Eo(t), Ei(t)) and (uo(t),ui(t)) denote the quasi- 
energies and Floquet states of the time-periodic Hamil- 
tonian 29] . Under the approximation we expect the fol- 
lowing relation, 

(7) 



where AEqi means quasi-energy splitting of the ground 
state doublet due to the tunneling between the wells. 

Let us confirm the relation in Eq.(7) numerically. In 
Fig. 5 we show the e— dependence of tv/AEqi- The be- 
havior is analogus to the e— dependence of the tunneling 
period of the oscillation Ph{t) in Fig. 4, in the weak 
perturbation regime. 

The similar correspondence between the tunneling pe- 
riod and the change of the quasi-energy splitting have 
been reported for the other double- well system by Tomso- 
vic et al. j23, El HS| • It is well-known that the chaos 
around the separatrix contributes to the enhancement 
of the tunneling split between the doublet, i.e. chaos- 
assisted tunneling. The reduction of the tunneling period 
can be approximately explained by the chaos-assisted in- 
stanton picture in the coherent oscilation regime e < e c . 
The simple explanation for the perturbed instanton pic- 
ture based on the width of the chaotic layer in the classi- 
cal dynamics is given in appendix B. (See also appendix 
A.) ' 
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FIG. 5: Perturbation strength dependence of quasi-energy 
splitting in the coherent motion regime e < e c in the 
monochromatically perturbed case (Af = 1). 



Generally speaking, as the number of frequencies M 
increases the tunneling period is more reduced as seen 
in Fig. 4 although we do not have analytic representation 
in the polychromatically perturbed cases. We conjecture 
that as seen in appendix A the increasing of the width 
of the stochastic layer contributes the reduction of the 
tunneling period even in the polychromatically perturbed 
cases. 
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FIG. 6: Plots of the uncertainty product AqAp as a function 
of time for some combinations of the parameters, (a) M — 2, 
e = 0.1. (b) M = 2, e = 1.0. (c) M = 10, e = 0.3. (d) 
M = 10, e = 1.0. 



C. uncertainly product 

Here let us investigate the spread of the wave packet in 
the phase space (q,p). Hitherto we mainly investigated 
the dynamics in q— space by Ph{t)- The phase space vol- 
ume gives a part of the compensating information for 
the phase space dynamics of the wave packet. Figure 6 
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FIG. 7: (a)Tunneling probability as a function of time 

t under stochastic perturbation with e = 0.4. (b)Plots of the 
uncertainty product AqAp versus time for various e's with the 
stochastic perturbation. The stochastic perturbation strength 
e is normalized to be equivalent to one of the polychromatic 
perturbation. 



presents the uncertainly product, i.e. phase space vol- 
ume, as a function of time for various cases, which is 
defined by, 



AqAp = V((q-(q)) 2 W((p-(j>))* 



(8) 



where (...) denotes quantum mechanical average. The un- 
certainty product can be used as a measure of quantum 
fluctuation [34| . The initial value is AqAp — h/2(= 0.5) 
for the Gaussian wave packet. It is found that in the case 
M = 2 the increase of the perturbation strength does not 
break the coherent oscillation and enhances the frequency 
of the time-dependence of the uncertainty product. For 
the relatively large e in M — 10, AqAp increases un- 
til the wave packet is relaxed in the space, and it can 
not return to Gaussian wave packet anymore. For the 
larger time scale, it fluctuates around the corresponding 
certain level. We can expect that the structure of the 
time dependence well corresponds to the behavior of the 
transition probability Pl{P) m Fig- 1- 



D. stochastically perturbed case 

It will be instructive to compare the above irregular 
motion under the polychromatic perturbation with the 
stochastically perturbed one. We recall that the stochas- 
tic perturbation, composed of the infinite number of the 
frequency components [M — > oo) with absolute continu- 
ous spectrum, can break the coherent dynamics. Indeed, 
if the time dependence of the potential comes up with 
the stochastic fluctuation as {{A{t\) — a)(A(t 2 ) — a)) en oc 
Ts5(t2 — t\), where (...) en and Tb denote ensemble av- 
erage and the temperature respectively, the stochastic 
perturbation partially models a heat bath coupled with 
the system 35] . Then the number of the frequency com- 
ponent corresponds to the number of degrees of free- 
dom coupled with the double- well system. The Pl (t) for 
the stochastic perturbation is shown in Fig. 6(a). The 
stochastic perturbation can be achieved numerically by 
replacing A(t) — a in the Eq.(2) by random number and 
we use uniform random number which is normalized so 
that the power of the perturbation is the same order to 
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one of the polychromatic case. In the limit of large M 
the motion under the polychromatic perturbation tends 
to approach the one driven by the stochastic perturbation 
provided with the same perturbation strengths e. 

Figure 7(b) shows the uncertainty product AqAp for 
the stochastically perturbed cases. It is found that 
the time-dependence of the uncertainty product in the 
stochastically perturbed case behaves similarly to the 
polychromatically perturbed ones for the relatively small 
e(<~ 0.1). On the other hand, for the relatively larger 
e(=0.4) the time-dependence shows quite different be- 
havior. While AqAp grows linearly with time in the 
stochastically perturbed case, in the polychromatically 
perturbed cases the growth of AqAp saturates at a cer- 
tain level. The linear growth of AqAp shows that the 
external stochasticity breaks the quantum interference in 
the internal dynamics. The growth of AqAp is strongly 
related to the growth of the energy of the packet [2(| . In 
the polychromatically perturbed cases the energy growth 
saturates at certain level due to quantum interference. 
On the other hand, in the case the energy grows un- 
boundedly, the activation transition becomes much more 
dominant than the tunneling transition when the wave 
packet transfers the opposite well. The details concern- 
ing relation between the stochastic resonance [35| and 
suppression of the energy growth will be given elsewhere 

Note that the polychromatic perturbation can be iden- 
tified with a white noise (or a colored noise if the frequen- 
cies are distributed over a finite band width) only in the 
limit of M — ► oo, while the stochastic perturbation can 
model a heat bath that breaks the quantum interference 
of the system. A similar phenomenon by the different 
property of the perturbation has been observed as "dy- 
namical localization" and the "noise-assisted mixing" of 
the quantum state in the momentum space in the quan- 
tum kicked rotor model fad . 



E. phase space portrait 

Finally, in order to see effect of quantum fluctuation, 
we compare the quantum states with the classical and 
semiquantal motions in the phase space for some cases. 
The semiquantal equation of motion is given by general- 
ized Hamilton-like equations as, 

dq 

dt = P > 
dp 

~dt 

d<t> 

~dt 

dm 

~dt 



s - (a) 

6 - 



A(t)q - 3h<j>q, 



4 - 
2 - 



(e) 



-4 - 
-6 - 




where the canonical conjugate pair (cj), w) is defined by 
the quantum fluctuation Aq and Ap as, Aq 2 = fuf), Ap 2 — 



FIG. 8: Classical ((a), (b)) and semiquantal ((c), (d)) stro- 
boscope phase space portrait of the polychromatically per- 
turbed double-well system with the same initial conditions 
(<lo,po) and the same parameters at t = 2im/f2i, n = 1, 2, .... 
The selection of initial condition of the fluctuation follows 
the minimum uncertainty. Contour plots ((e), (f)) of the 
Hushimi functions for the corresponding the quantum state 
at t — 8.7 x 10 3 . Contour lines in the panel (e) at the values 
0.01, 0.02, 0.05, 0.08 and 0.1, in the panel (f) at 0.01, 0.02, 
0.03, 0.04 and 0.05. M = 7 and e = 0.4 for (a), (c) and (e). 
M = 7 and e = 0.8 for (b), (d) and (f). 



ft(jr +47r 2 </>). For more details consult It is directly 
observed that the quantum tunneling phenomenon en- 
hances chaotic motion in comparing to the classical and 
semiquantal trajectories. In Fig. 8(a) and (b) Poincare 
surface of section of the classical trajectories in the phase 
plane at M — 7 are shown. The stroboscopic plots are 
taken at t = 2irn/ili(n — 1,2,...), due to non time- 
periodic structure of the Hamiltonian. In the relatively 
small perturbation strength e = 0.4, the trajectories stay 
the single well, and are stable even for the long-time evo- 
lution. Figure 7(c) and (d) show the Poincare section of 
the semiquantal trajectories for a polychromatically per- 
turbed double-well system with M = 7 (the stroboscopic 
plots are taken at t = 2~Kn/£l\ again). The semiquantal 
trajectories for the squeezed quantum coherent state can 
be obtained by an effective action which includes par- 
tial quantum fluctuation to all order in h 3]. It can be 
seen that in comparing with ones of classical dynamics 
the trajectories in the semiquantal dynamics spreads into 
the opposite well even for the small e. This corresponds 
to the quantum tunneling phenomenon through the semi- 
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quantal dynamics. Apparently, the partial quantum fluc- 
tuation in the semiquantal approximation enhances the 
the chaotic behavior. Notice that the semiquantal picture 
breaks down for the irregular quantum states because the 
centroid motion becomes irrelevant. 

In Fig. 8(e) and (f) the corresponding coherent state 
representation for the quantum states are shown. It is 
directly seen that the wave packet spreads over the two- 
wells and the shape is not symmetric. Once the wave 
packet incoherently spreads over the space, it can not 
return to the initial state anymore. We have confirmed 
that in a case without separatrix (single- well), namely 
the case that a in Eq. (2) is replaced by —a, in the clas- 
sical phase space the coherent oscillations have remained 
against the relatively large M and/or e. It follows that 
the full quantum interference suppresses the chaotic be- 
havior as seen in the semiquantual trajectories. 



IV. SUMMARY AND DISCUSSION 

We numerically investigated influence of a polychro- 
matic perturbation on wave packet dynamics in one- 
dimensional double-well potential. The calculated phys- 
ical quantities are the transition rate the time- 
fluctuation APl, uncertainty product AqAp and phase 
space portrait. The results we obtained in the present 
investigation are summarized as follows. 

(1) We classified the motions in the parameter space 
spanned by the amplitude and the number of frequency 
components of the oscillatory driving force, i.e. coherent 
motions and irregular motions. The critical value e c (M) 
which divides the behavior of Ph{i) into regular and ir- 
regular motions depends on the number of the frequency 
component M. 

(2) Within the regular motion range, the period of 
the tunneling oscillation is reduced with increase of the 
number of colors and/or strength of the perturbation. It 
could be explained by the increase of the instanton tun- 
neling rate due to appearance of the stochastic layer near 
separatrix |28| . In this parameter regime the perturbed 
instanton pict ure is one of expression for chaos-assisted 
tunneling |32j and chaos-assisted ionization picture re- 
ported for some quantum chaos systems |33| . 

(3) In the irregular motion in the polychromatically 
perturbed cases, the growth of AqAp initially increases 
and saturates at certain level due to quantum interfer- 
ence. On the other hand, in the stochastically perturbed 
case the uncertainty product grows unboundedly because 
the external stochasticity breaks the quantum interfer- 
ence in the internal dynamics. The growth of AqAp is 
strongly related to the growth of the energy of the wave 
packet. 

(4) It is expected that the quantum fluctuation are 
always large for the classically chaotic trajectories com- 
pared to the regular ones. This implies that the quantum 
corrections to the evolution of the phase space fluctua- 
tion become more dominant for classically chaotic trajec- 



tories. 

(5) In the semiquantal approximation the partial quan- 
tum fluctuation enhances the chaotic behavior, and si- 
multaneously the chaos enhances the tunneling and de- 
coherence of the wave packet. The quantum fluctuation 
observed in the semiquantal picture is suppressed by in- 
terference effect in the fully quantum motion. The semi- 
quantal picture can not apply to the chaos-induced delo- 
calized states. 

Furthermore, in the appendices, we gave classical 
phase space portraits in the polychromatically perturbed 
double-well system and a simple explanation for the per- 
turbed instanton tunneling picture for the reduction of 
the tunneling period in the coherent motion regime. 



Although we have dealt with quantum dynamics of 
wave packet with paying attention to existence of the en- 
ergetic barrier, we can expect that the similar phenomena 
would appear by dynamical barrier in the system. The 
details will be given elsewhere |26j |. 

APPENDIX A: EFFECT OF POLYCHROMATIC 
PERTURBATION ON CLASSICAL PHASE 
SPACE PORTRAITS 

We show classical stroboscopic phase space portrait 
in this appendix with paying an attention to the effect 
of polychromatic perturbation on the chaotic behavior. 
In the classical dynamics, such a system shows chaotic 
behavior by the oscillatory force A(t) 0i • The New- 
ton's equation of the motion is 

fl-A(t)q + q 3 = 0. (Al) 

Note that in the monochromatically perturbed case 
(M = 1, A(t) = a — e sin fit), the equation is known 
as nonlinear Mathieu equation which can be derived 
from surface acoustic wave in piezoelectric solid |37J| and 
nanomechanical amplifier in micronscale devices [381 ]. 



1. Classical phase space potraits 

In Fig. A.l, we show the change of the classical stro- 
boscopic phase space portrait changing the perturbation 
parameters. Increasing the perturbation strength e de- 
stroys the separatrix and forms a chaotic layer in the 
vicinity of the separatrix. Needless to say, the phenom- 
ena have been observed even in the monochromatically 
perturbed cases 24] . In the polychromatically perturbed 
cases {M > 1) the smaller the strength e can generate 
chaotic behavior of the classical trajectories the larger 
M is j2(J . It should be emphasized that in the polychro- 
matically perturbed cases the width of the chaotic layer 
grows faster than the monochromatically perturbed case 
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as the perturbation strength increases. As a result, the 
increase of the color contributes the increase of the width 
of the stochastic layer in the polychromatically perturbed 
cases. 
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FIG. 10: Maximal Lyapunov exponents \™ ax estimated by 
some classical trajectories within the finite-time t € [0, T], 
where T is tunneling time given inSubsect.3.1. 



We can roughly observe a transition from motion of 
KAM system (\™ l ax ~ 0) to chaotic motion (\% ax > 0) 
as the perturbation strength increases. The increasing of 
the number of color M reduces the value of the critical 
perturbation strength e~ of a transition from a motion 
of KAM system to fully chaotic motion. Roughly speak- 
ing, a transition of the classical dynamics corresponds to 
the transition from coherent motion to irregular one in 
the quantum dynamics. We expect that the transition 
observed in Sect. 3 will corresponds to quantum signa- 
tures of the KAM transition from the regular to chaotic 
dynamics. 



Some considerations to the effect of 
polychromatic perturbation 



FIG. 9: Stroboscopic classical phase space portraits of the 
polychromatically perturbed double-well system. Cross sec- 
tions of (p, q) space are plotted at 2irn/Qi,n = 1, 2, 3, .... The 
e— dependence for M = 5 is shown in (a)e = 0.1, (b)e = 0.3, 
(c)e = 0.4 and (d)e = 0.8. The M-dependence for e = 0.5 is 
shown in (e)M = 1, (f)M = 2, (g)M = 5 and (h)M = 8. 

Here, we use the increasing rate of infinitesimal dis- 
placement along the classical trajectory for the extent 
of chaotic behavior as a finite-time Lyapunov exponent 
\™ ax . We prepare various initial points in the phase 
space, and conveniently adapt a trajectory with max- 
imal increasing rate among the ensemble within the 
finite-time interval as the finite-time Lyapunov exponent. 
Note that an exact Lyapunov exponent should be de- 
fined for the long-time limit. However, the roughly esti- 
mated Lyapunov exponent is also useful to observe the 
classical-quantum correspondence. Figure A. 2 shows the 
e— dependence of classical Lyapunov exponents for vari- 
ous cases estimated by the numerical data of the classical 
trajectories. 



In this subsection, we conceptually consider a role of 
the polychromatic perturbation different from monochro- 
matic one. Note that the change of the number M of 
colors also changes the qualitative nature of the under- 
lying dynamics because M corresponds to the effective 
number of DOF under some conditions El). In the 
our model, when the number of DOF of the total system 
is more than four, i.e. M > 3, the classical trajectories 
can diffuse along the stochastic layers of many resonances 
that cover the whole phase space if the trajectory starts 
in the vicinity of a nonlinear resonance. The number 
of resonances increases rapidly with DOF, changing the 
characteristic of population transfer from bounded to dif- 
fusive. Such a global instability is known as Arnold dif- 
fusion in nonlinear Hamiltonian system with many DOF 
|39|| . The effect of Arnold diffusion in quantum system is 
not trivial and the study just has started recently |40| . 
The more detail is out of scope of this paper. 

Moreover, we can regard the time-dependent model 
of Eq.(l) as nonautonomous approximation for an au- 
tonomous model, consisting of the double-well system 
coupled finite number M of harmonic oscillators with 
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the incommensurate frequencies {Qj,j — 1,1,. ..,M}. It 
is worth noting that the linear oscillators can be iden- 
tified with a highly excited quantum harmonic oscilla- 
tors, which all phonon modes are excited around Fock 
states with large quantum numbers. Then the above 
model can be regarded as a double-well system coupled 
with M phonon modes. Without the interaction with the 
phonon modes the Gaussian wavepacket remains the co- 
herent motion. Then the number of DOF of total system 
is 1 + M and the number of the frequency components 
M corresponds to the that of the highly excited quantum 
harmonic oscillators. The detail of the correspondence is 
given in Ref. p5| . 

In quantum chaotic system with finite and many DOF 
we expect occurrence of a dissipative behavior. For exam- 
ple, we consider simulated light absorption by coupling a 
system in the ground state with radiation field. Then sta- 
tionary one-way energy transport from photon source to 
the system can be interpreted as occurrence of quantum 
irreversibility in the total system. Such a irreversibility is 
called chaos-induced dissipation in quantum system with 
more than two DOF [4l|. In this sense, we can expect 
occurrence of the one-way transport phenomenon in the 
delocalized state in the irregular motion phase if it cou- 
ples with the other DOF in the ground state as seen in 
Ref. HI. 



1. A discrete mapping 

Indeed, trajectories in the neighborhood of the sepa- 
ratrix of the system are well reproduced by the whisker 
map of the system. Whisker map is a map of the en- 
ergy change I n and phase change </>„ of a trajectory in 
the neighborhood of the separatrix for each of its motion 
during one period of the perturbation, i.e. action-angle 
variable. Moreover, if we linearize the whisker map which 
describes the behaviors of the trajectories in the neigh- 
borhood of the fixed point, we can obtain the following 
standard map, 

f I n +i = I n - K(e, fi) sin<^„ , , 

\ <t>n+l = <t>n + I n +1 

, where K (e, fi) is a nonlinear parameter of local instabil- 
ity that the exact function form which is not essential for 
our purpose. K increases with e, and K > 1 means that 
the dynamics of the system is locally unstable. A com- 
parison has done between the whisker map and the strobe 
plots in the time continuous version by Yamaguchi [2l| . 
The form of the mapping is convenient for the estimate 
of the width of the stochastic layer. 

Kuvshivov et al has given the estimated parameter of 
local instability, width of stochastic layer, and correlator 
for perturbed instanton solution for the perturbed pen- 
dulum system |2S| . 



APPENDIX B: PERTURBED INSTANTON 
TUNNELING 

In this appendix, we consider the reduction of the tun- 
neling period as the perturbation strength increases in 
the coherent oscilation regime e < e c , based on a per- 
turbed instanton tunneling. In a double- well system with 
dipole-type interaction, qsm(Qt), the energitical barrier 
tunneling between the symmetric double-well can be ex- 
plained by a three-state model or chaos-assisted tunnel- 
ing (CAT) US IH H H3- The three states that take 
part in the tunneling are a doublet of quasi-degenerate 
states with opposite parity, localized in the each well, 
and a third state localized in the chaotic layer around 
the separatrix. 

However, note that less attention has been paid to tun- 
neling in KAM system while chaotic dynamics has been 
modeled by multi-level Hamiltonian and random ma- 
trix model to describe the chaos-assisted tunneling. We 
give an expresion of the tunneling amplitude in " chaos- 
assisted instanton tunneling" firstly proposed by Kuvshi- 
nov et al for a Hamiltonian system with time-periodic 
perturbation |2g. Let us consider only monochromat- 
ically perturbed case (M = 1 in Eq.(2), Q, — fii) be- 
cause the separatrix destruction mechanism by the time- 
periodic perturbation has an universality although our 
system is different from their one. . 



2. Action of the instanton tunneling 

The perturbation destroys separatrix of the unper- 
turbed system and the stochastic layer appears. In the 
regular motion denoted by the circles in Fig. 3, classical 
chaos can increase the rate of instanton tunneling due 
to appearance of the stochastic layer near separatrix of 
the unperturbed system. As a result the frequency of 
time-dependence Pz,(t) increases as the classical chaos 
becomes remarkable in the parameter regime. Note that 
the perturbed instanton tunneling picture disappears in 
the strongly perturbed regime due to the derealization 
of wavepacket. 

Here we give only relation between the width of the 
stochastic layer and the tunneling amplitude in terms of 
path integral in imaginary time r, found by Kuvshivov 
et al. Tunneling amplitude between the two wells in the 
perturbed system can be given by integration over energy 
of tunneling amplitude A tU n in unperturbed system as, 

A tun = / dE D[q(T,E)]exp{-S[q(T,m^ 

J0 J q(r)=— g 

where S[q(r, E)\ denotes the Euclidian action. AH = 
2\H S — Hb\ denotes the width of stochastic layer, where 
H s and Hf, are the energy of the unperturbed system on 
the separatrix and on the bound of stochastic layer, re- 
spectively. q(r, E) is classical solution of Euclidian equa- 
tion of motion. The contribution of the chaotic instanton 
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solution are taken into account by means of integration 
over E which is energy of the instanton. 

The perturbed instanton solutions correspond to the 
motions in vicinity of the separatrix inside the layer. 
The only manifestation of the perturbation in this ap- 
proximation is the appearance of a number of additional 
solutions of the Euclidian equation of motion with energy 
close to the energy of the unperturbed one-instanton so- 
lution inside the stochastic layer. Accordingly we can 
expect that the appearance of the stochastic layer en- 
hances the tunneling rate as reported in the other sys- 
tems |28| . However, we have to have in mind that the 
result is obtained in the first order on coupling constant 
e of the time-periodic perturbation and does not take 
into account the structure of stochastic layer. The ap- 
proximation is valid if the layer is narrow by neglecting 
the higher order resonances in the phase space. For the 
more details of the perturbed instanton see Ref . |2^| . The 



increasing of the tunneling amplitude is directly related 
to energy splitting A_B i between the quasi-degenerate 
ground Floquet states. 

As seen in Fig.A.l, the increasing of number of color 
M can enhance the width of the stochastic layer with 
the perturbation strength being kept at a constant value. 
The theoretical explanation for the reduction of the tun- 
neling period with the number of color is open for further 
study. 

We expect that in the double-well system under the 
polychromatic perturbation this numerical study will be 
useful for the analytical derivation of "reduction of tun- 
neling period" and " critical strength of a transition from 
localized to delocalized behavior of wavepacket" by ex- 
tension of the monochromatically perturbed case. The 
chaos assisted instanton theory might be applicable if we 
will exactly estimate the width of the stochastic layer in 
the system under the polychromatic perturbation. 
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